Crossing the dripline to n N using elastic resonance scattering 
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The level structure of the unbound nucleus n N has been studied by 10 C+p elastic res- 
onance scattering in inverse geometry with the LISE3 spectrometer at GANIL, using 
a C beam with an energy of 9.0 MeV/u. An additional measurement was done at 
the A1200 spectrometer at MSU. The excitation function above the 10 C+p threshold 
has been determined up to 5 MeV. A potential-model analysis revealed three resonance 
states at energies 1.27±S:S| MeV (F=1.44±0.2 MeV), 2.0llo;os MeV, (r=0.84±0.2 MeV) 
and 3.75±0.05 MeV, (F=0.60± 0.05 MeV) with the spin-parity assignments P=l/2 + , 
l/2~, 5/2 + , respectively. Hence, n N is shown to have a ground state parity inversion 
completely analogous to its mirror partner, n Be. A narrow resonance in the excitation 
function at 4.33±0.05 MeV was also observed and assigned spin-parity 3/2~. 

PACS number(s): 21.10.Hw, 21.10.Pc, 25.40. Ny, 27.20. +n 



I. INTRODUCTION 

The exploration of exotic nuclei is one of the 
most intriguing and fastest expanding fields in 
modern nuclear physics. The research in this do- 
main has introduced many new and unexpected 
phenomena of which a few examples are halo sys- 
tems, intruder states, soft excitation modes and 
rare /3-delayed particle decays. To comprehend 
the new features of the nuclear world that are re- 
vealed as the drip-lines are approached, reliable 
and unambiguous experimental data are needed. 
Presently available data for nuclei close to the 
driplines mainly give ground-state properties as 
masses, ground state J 7r and beta-decay half-lives. 
Also information on energies, widths and quantum 
numbers I* of excited nuclear levels are vital for 
an understanding of the exotic nuclei but are to a 
large extent limited to what can be extracted from 



(3 decays. Nuclear reactions can give additional 
information, in particular concerning unbound nu- 
clear systems. However, the exotic species are 
mainly produced in complicated reactions between 
stable nuclei. These processes are normally far too 
complex to allow for spin-parity assignments of the 
populated states, and hence are of limited use for 
spectroscopic investigations. Instead of using com- 
plex reactions between stable nuclei, the driplines 
can be approached in simple reactions involving ra- 
dioactive nuclei. An example is given in this paper 
where elastic resonance scattering of a 10 C beam 
on a hydrogen target was used to study the un- 
bound nucleus n N. With heavy ions as beam and 
light particles as target, the technique employed 
here is performed in inverse geometry. The use 
of a thick gas target instead of a solid target is 
another novel approach. This technique has been 
developed at the Kurchatov Institute H where it 
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has been employed to study unbound cluster states 
with stable beams [||. The perspectives of using 
radioactive beams in inverse kinematics reactions 
to study exotic nuclei are discussed in and the 
method was used in ||] . Resonance elastic scatter- 
ing in inverse kinematics using radioactive beams 
and a solid target has been used at Louvain-la- 
Neuve §§. 

This experiment is part of a large program for 
investigating the properties of halo states in nu- 
clei (jfj. A well studied halo nucleus is 11 Be where 
experiments have demonstrated that the ground 
state halo mainly consists of an lsi/2 neutron cou- 
pled to the deformed 10 Be core [|]||, in contra- 
diction to shell-model which predicts that the odd 
neutron should be in a Opi/2 state. The Opi/ 2 level 
is in reality the first excited state, while the ground 
state is a lsi/2 intruder level JhJ. This discov- 
ery has been followed by numerous papers investi- 
gating the inversion, e. g. references 0@. The 
mirror nucleus of 11 Be, n N, should have a 1/2+ 
ground state with the odd proton being mainly in 
the lsi/2 orbit, if the symmetry of mirror pairs 
holds. However, n N is unbound with respect to 
proton emission which means that all states are 
resonances that can be studied in elastic scatter- 
ing reactions. The first experiment devoted to a 
study of the properties of the low-energy structure 
of 11 N used the three-nucleon transfer-reaction 
14 N( 3 He, 6 He) n N. The results indicated a resonant 
state at 2.24 MeV which was interpreted as 
the first excited l/2~ state rather than the 1/2+ 
ground state. 

In this paper we present excitation functions 
at laboratory (lab.) angles of 0° measured at 
GANIL |J and MSU, and at 12.5° with respect to 
the beam direction measured at GANIL. A thor- 
ough analysis, using a potential model as well as 
a simplified R-matrix treatment, gives unambigu- 
ous determination of the quantum structure of the 
three lowest resonances in the 10 C+p system. 



limits the applicability of this method to investi- 
gations in regions close to (3 stability. However, 
it is possible to produce dripline species in simple 
reactions involving radioactive nuclei. When using 
this approach, the beam is composed of radioac- 
tive ions and the target of light nuclei, eliminating 
the need for a radioactive target. Since this is the 
inverse setup to the one traditionally used in scat- 
tering experiments, the method is usually denoted 
as elastic scattering in inverse geometry. 

The advantage of using gas instead of a solid 
target is twofold. Firstly, the thickness of a gas 
target can be changed continuously and easily by 
adjusting the gas pressure, and secondly the tar- 
get is very homogeneous. The beam parameters 
of radioactive ion beams (RIB's) are limited; the 
spread in both energy and space are much larger 
than what can be obtained for stable beams, and 
the intensities are of course much smaller. As 
will be seen below, the beam properties are not 
of great importance in the experimental approach 
used here. Elastic resonance scattering is charac- 
terized by large cross sections and is therefore well 
suited for use with low-intensity RIB's. These and 
several other features of the elastic resonance scat- 
tering in inverse geometry on thick targets will be 
illuminated in the following subsections. 

The expressions for elastic cross section in the 
case of proton scattering on spin-less nuclei, eqs. [l] 
and |^, can for example be found in |l6| , 
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II. ELASTIC RESONANCE SCATTERING: 
METHOD AND FORMALISM 

The first description of elastic resonance scat- 
tering was given by Breit and Wigner and 
it is now a theoretically well understood reaction 
mechanism [jlqJlq] . Traditionally, elastic resonance 
scattering experiments have been performed by 
bombarding a thin target with a light ion beam, 
narrow in space and time. To obtain an excitation 
function the beam energy then had to be changed 
in small steps of the order of the experimental res- 
olution. The need for a radioactive target severely 
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Symbols introduced here are defined as: 
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± : denotes states with j = I ± | 

z: charge of the proton 

Z: charge of the spin-zero particle 

/i: reduced mass 

v: the relative velocity of the particles 

k: the magnitude of the wave vector 

at- the Coulomb phase shift 

Pl(cos6): Legendre polynomial 

P} {cosff) : associated Legendre polynomial 

The first term in A(9) represents the Coulomb 
scattering. The other terms in A{9) and B(9) ex- 
press scattering due to nuclear forces. The phase 
shift 81 is the sum of the phase shift from hard 
sphere scattering, —<j>i % and the resonant nuclear 
phase shift, (3e,: 
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The differential cross-section has its maximum 
in the vicinity of the position where the phase shift 
passes through (n + |)7r. Therefore, a frequently 
used definition of the resonance energy is where 6 — 
(n + i)7r, see section IV. It is favourable to study 
resonance scattering at 180° cm. where eq. [j] is 
simplified. At this angle, only m = sub-states 
contribute to the cross section and both potential 
and Coulomb scattering are minimal. An advan- 
tage of the inverse geometry setup is its possibility 
to measure at 180° cm. 



A. Kinematical relations 

We define the laboratory energies of the bom- 
barding particles before the interaction in inverse 
(E) and conventional (T) geometry as Eq and To, 
respectively. The notation used mainly follows 
ref. |rH| , primed energies being in the cm. sys- 
tem: 

m,M: mass of the light and heavy particles 
Em, Tm'- heavy particle M laboratory energies 

after interaction 
E m , T m : light particle m laboratory energies 

after interaction 
Oiab. ■ scattering angle of the light particle in 

the laboratory system 



The relations between laboratory energy of the 
beam and the cm. energy of the heavy nucleus: 
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The expressions for the lab. energies of the light 
particle m that will be detected after scattering: 
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In the equation above, K is the ratio of the masses 
{E /T =M/m = K since E M = T M ). Inserting 
Oiab. = 0° in eq. || leads to the following ratio be- 
tween the energy of the measured particle in con- 
ventional and inverse geometry: 
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As is seen from eq. [?], the detected energy of the 
light particles is close to 4 times higher for inverse 
kinematics as compared to the conventional geom- 
etry at the same cm. energy. This is an impor- 
tant gain for the study of resonant states near the 
threshold. The excitation energy in the M + m 
compound system is obtained as the sum of the 
cm. energies for particles m and M: 
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Using eq. ^, this can be expressed in terms of the 
measured particle energy E rn . In case of inverse 
kinematics, the excitation energy of the compound 
system becomes: 
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Because of the low energies involved, a non- 
relativistic expression can be used. 



B. General set-up of elastic scattering in 
inverse geometry 

The basic experimental setup consists of a ra- 
dioactive ion beam which is incident on a scatter- 
ing chamber filled with gas. The thickness of the 
gas target is adjusted to be slightly greater than 
the range of the beam ions. Charged-particle de- 
tectors are placed at and around the beam direc- 
tion, i.e. 180° cm., as shown in Figure [l]. As 
they are continuously slowed down in the gas, the 
beam ions effectively scan the energy region from 
the beam energy down to zero, giving a continuous 
excitation function in this interval. When the en- 
ergy of the heavy ion corresponds to a resonance in 
the compound system, the cross section for elastic 
scattering increases dramatically and can exceed 
1 b, making it possible to neglect the non-resonant 
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contributions which are on the order of nib. For the 
ideal case of a mono-energetic beam, each interac- 
tion point along the beam direction in the chamber 
corresponds uniquely to one resonance energy and, 
as we study elastic scattering, to a specific proton 
energy for each given angle. Because the distance 
from the detector is different for each proton en- 
ergy, the solid angle also varies with proton energy 
and is quite different for low and high energy res- 
onances. 

The high efficiency of the method is mainly a 
result of the large investigated energy region. If 
we compare the scanned region of 5-10 MeV with 
the typical energy step of wlO-20 keV in conven- 
tional scattering measurements, the gain is 250- 
1000 times. 



C. Energy resolution 

The initial energy spread of our 10 C beam was 
1.5% of the total energy, which naturally increased 
along the beam path in the gas. The energy spread 
of the beam results in excitation of the same reso- 
nance at different distances from the detector. As- 
suming that AE is the energy spread at some point 
in the gas, this distance interval Aa; is given by 
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where (4^)m is the specific energy loss of the beam 
nuclei in the gas. Due to the protons energy loss in 
the gas, the measured proton energies correspond- 
ing to the same resonance are slightly different. 
The resulting spread of proton energies, e, corre- 
sponding to the interval Ax will be 



e = AE 



V dx > m 
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Here, (4^) m denotes the specific energy loss of the 
recoil nuclei (protons) in the gas. Taking into ac- 
count the different velocities of the beam ions and 
the scattered protons as well as the Bethe-Bloch 
expression for specific energy loss, one finds: 



AE z 2 
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In the case of 10 C+p interaction, eq. [l2] becomes 
Hence, for AE = 5 MeV a lab. energy reso- 
lution of 35 keV is expected. The effective cm. en- 
ergy resolution will be about four times better than 
the resolution in the lab. frame, see eq. 0. Thus it 
is clearly shown that the energy spread of the ra- 
dioactive beam does not restrict the applicability 
of the method. Many other factors influence the 



final resolution, for example the size of the beam 
spot and detectors, the detector resolution, the an- 
gular divergence of the beam and straggling of light 
particles in the gas. These factors can be taken into 
account by Monte Carlo simulations. In reality, an 
effective energy resolution of 20 keV in the cm. 
frame is feasible. At angles other than 180° the 
resolution deteriorates, mainly due to kinematical 
broadening of the energy signals for protons scat- 
tered at different angles. This contribution to the 
resolution could be reduced by tracking the proton 
angles. 



D. Background sources 

A cornerstone of the described experimental ap- 
proach is that elastic resonance scattering domi- 
nates over other possible processes. The compet- 
ing reaction channel which has to be treated for 
each specific case is inelastic resonance scattering, 
as it is a resonant process which produces the same 
recoil particles as the elastic scattering. However, 
the elastic and inelastic resonance scattering reac- 
tions can be distinguished from each other. The 
energy of the scattered nuclei from inelastic reso- 
nance scattering at 0° is given by eq. p~3] if E* / E <g^ 
1, where E* is the excitation energy of the beam 
nucleus Jl7[ . 
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Comparing this with eq. |^, one sees that the en- 
ergy of heavy ions has to be larger by an amount £ 
for the inelastic scattering to obtain the same en- 
ergy of a light recoil from the elastic and inelastic 
scattering reactions, when £ is defined in eq. |l4| . 
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For the 10 C+p case, where £*( 10 C(2+)) = 
3.35 MeV, eq. |l4| shows that the inelastic resonance 
scattering should take place at about 20 MeV 
higher energy than the elastic one for the two pro- 
cesses to mix in the elastic scattering excitation 
function. The inelastic resonance reaction thus has 
to take place further from the detectors, closer to 
the entrance window, in order to produce a scat- 
tered particle with the same energy as the cor- 
responding elastic process. The two processes in 
question hence can give the same energies of the re- 
coil protons but their ToF (window-detector) will 
differ. The time difference between the two types 
of events will be on the order of a few ns, and can 
thus be separated in the analysis. No such events 
were seen in our data. 
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Other scattering reactions contribute very little 
to the spectrum, especially at 180° cm., the excep- 
tion being low energies where the Coulomb scatter- 
ing cross sections increases. However, this scatter- 
ing is well understood and can be included in the 
data treatment. Additional sources of background 
are (3 particles from decaying radioactive ions in 
the gas, beam ions which penetrate the gas target, 
and particles scattered in the entrance window. 

III. EXPERIMENTAL PROCEDURE 

The first experiment was performed using the 
LISE3 spectrometer at the GANIL heavy-ion fa- 
cility. The secondary 10 C beam was produced 
by a 75 MeV/u 12 C 6+ beam with an intensity of 
2-10 12 ions/s which bombarded a 8 mm thick, ro- 
tating Be target and a fixed 400 jiim Ta target. 
The 10 C fragments were selected in the LISE3 
spectrometer, using an achromatic degrader at 
the intermediate focal plane (Be, 220 /im thick) 
and the Wien-filter after the last dipole. The 
50 cm long scattering chamber was placed at the 
final focal plane. Immediately before the 80 /im 
thick kapton entrance window, a PPAC (Parallel 
Plate Avalanche Counter) registered the position 
of the incoming ions. The intensity of the sec- 
ondary beam, measured by the PPAC, was approx- 
imately 7000 ions/s, and due to the degrader and 
Wien filter a very low degree of contamination was 
achieved. The efficiency of the PPAC at this inten- 
sity and ion charge is close to 100%, which makes 
it easy to use the PPAC count rate to obtain ab- 
solute cross-sections. The scattering chamber was 
filled with CH4 gas, acting as a thick proton tar- 
get for the incoming 10 C ions. The gas pressure 
was adjusted to 816±5 mbar, which was the pres- 
sure required to stop the incoming beam just in 
front of the central detector. It is desirable to 
stop the beam close to the detectors in order to 
avoid loosing any protons scattered from a possi- 
ble low-lying resonance in n N. In the far end of the 
chamber an array of Si-detectors was placed. The 
detectors had diameters of 20 mm and thickness 
of 2.50 mm, corresponding to the range of 20 MeV 
protons. The time between the radio frequency 
(RF) from the cyclotron and the PPAC gave one 
time-of- flight signal (ToF), while the time differ- 
ence of the PPAC and detector signal gave ad- 
ditional ToF-information. The complete setup is 
shown in Figure [I]. 

As a first measurement, a low intensity 10 C 
beam was sent into the evacuated scattering cham- 
ber to get the total energy and spread of the sec- 
ondary beam after the foil, and this was deter- 



mined to be 90 MeV with a FWHM = 1.5 MeV. For 
background measurements, the scattering cham- 
ber was filled with CO2 gas at 450±5 mbar and 
bombarded with 12 C and 10 C beams, respectively. 
For our purposes, we assume that 16 and 12 C be- 
have similarly in proton scattering reactions. The 
measurements with the CO2 target would reveal 
any background stemming from the carbon nuclei 
in the CH4 target gas or from the kapton win- 
dow. Beam contaminations would also be present 
in these runs, and those background sources can 
subsequently be subtracted from the experimental 
excitation functions. 




a source 

FIG. 1. A schematic view of the scattering chamber 
and Si-detectors. The 12.5° detector has this angle to 
the middle of the chamber. The interaction position 
and the corresponding angle is calculated for each event 
when analyzing the data. 

The standard beam diagnostics observed admix- 
tures of d, a and 6 Li with the same velocity as the 
12 C secondary beam, while no contaminant parti- 
cles could be seen in the 10 C beam. The 10 C+CO2 
spectrum showed no prominent structure and was 
found to contribute less than 10% to the total cross 
section. This background spectrum was subtracted 
from the 10 C+CH4 spectrum before transforma- 
tion to the cm. system. 

Since 10 C is a f3 + emitter with a half life of 19.3 s, 
it is necessary to discriminate the positron sig- 
nals from the protons. This was done by selecting 
the protons in a two-dimensional spectrum show- 
ing ToF (PPAC to Si-detector) versus detected en- 
ergy, where the positrons are clearly distinguished 
from protons both by their uniform time distribu- 
tion and their maximum energy of 1.93 MeV. A 
positron with energy in this interval has a maxi- 
mum energy loss of 1.25 MeV in 2.50 mm Si, which 
simulates a scattered proton energy of 0.344 MeV 
in the excitation function of 11 N. Since the positron 
energies are small enough to lie in the energy range 
of Coulomb scattered protons, cutting away all 
events below this energy does not distort the in- 
teresting parts of the proton spectrum, as is seen 
in the inset in Figure |3| 

In this paragraph we justify our ignoring the 
background contributions to our spectra from in- 
elastic scattering of 10 C on hydrogen with excita- 
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tion of the particle stable 2 + level at 3.35 MeV 
in 10 C. The contribution from inelastic scattering 
has been estimated using available data on inelas- 
tic scattering of protons on a 10 Be target [JOl and 
a DWBA extrapolation to the whole investigated 
interval of energies. This shows that the contribu- 
tion from inelastic scattering does not exceed 1% 
of the observed cross section. 

The energy calibration of the Si detectors was 
done with a triple a-source ( 244 Cm, 241 Am and 
239 Pu) which was placed on a movable arm inside 
the chamber. Another calibration, and at the same 
time a performance check of the setup, was done 
by investigating known resonances in 13 N. The pri- 
mary 12 C beam, degraded to 6.25 MeV/u, was 
scattered on the methane target using a gas pres- 
sure of 240±5 mbar. The resulting proton spec- 
trum clearly shows the two closely lying resonances 
in 13 N (3.50 MeV, width 62 keV, and 3.55 MeV, 
width 47 keV, as can be seen in Figure H. 
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FIG. 2. Experimental spectrum of protons from the 
scattering of 12 C. The energy is given as c.m. energy 
of 13 N. The peaks are fitted with two coherently added 
resonances as described in the text. 

These resonances are overlapping and the width 
of the peak is 50 keV. The solid curve in Figure || 
is a fit obtained by coherently adding two curves in 
order to take interference into account. The 5/2 + 
resonance at 3.55 MeV has single particle (SP) 
nature Jl9| and was described using the potential 
model outlined in section while a Breit-Wigner 
curve was used for the 3/2~ state at 3.50 MeV. 
The resonant 1/2+ state in 13 N, 420 keV above 
the 12 C+p threshold, is not seen as it is overlap- 
ping the Coulomb scattering which dominates be- 
low 0.5 MeV. From the calibration measurements 
described above an energy resolution of 100 keV 
in the lab. frame was deduced, mainly determined 
by the detector resolutions and proton straggling 
in the gas. 

The experimental proton spectrum was, after 
subtraction of the measured carbon background, 



transformed into differential cross-section as a 
function of the excitation energy of n N, in the fol- 
lowing referred to as the excitation function of 11 N. 
Since each interaction point along the beam di- 
rection ideally corresponds to a specific resonance 
energy, the measured proton energy can, after cor- 
rection for its energy loss in the gas, be used to 
find the resonance energy in n N in the c.m. sys- 
tem. The inset in Figure ^| shows the experimental 
data as measured proton energy versus counts at 0° 
before the corrections for solid angle was made. 
Comparing this pictures to the one obtained after 
transformation to the c.m. system clearly shows 
the effect of differing solid angles for different pro- 
ton energies. The cross sections in the high energy 
part increases relative to the low energy part, as 
is clearly seen when comparing the inset to the 
transformed spectrum in Figure ||. 

Extracting the cross section from the data is 
straightforward, and the transformation to c.m. is 
done using eq. [Iq. 
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The relation between the scattering angle in the 
lab. and c.m. systems is simply 6> c . m . = 180°-2#; a fc.. 
The excitation function obtained after background 
subtraction and transformation into the c.m. sys- 
tem is shown in Figure ^|. 
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FIG. 3. The excitation function of n N is shown af- 
ter background subtraction and transformation to the 
c.m. system. The filled black circles represent the 
GANIL data, and the white squares show the result of 
the MSU experiment. The inset shows the raw data 
in the lab. system. The upper scale in the inset is 
a rough calibration to c.m. energy, shown for com- 
parison, while the lower scale is the detected proton 
energy. 

The more detailed analysis now performed re- 
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vises the absolute cross section to a larger value 
from what was previously published in M . 

A second independent measurement of 10 C+p 
elastic scattering using the same method was made 
at NSCL where the A1200 spectrometer delivered 
the 10 C beam. The experimental conditions were 
the same as in the GANIL experiment, except that 
at NSCL a AE-E telescope was placed at 0° and no 
Wien-fhter was used. The energy of the 10 C beam 
after the foil was 7.4 MeV/u and the beam inten- 
sity was 2000 pps. The data from these two ex- 
periments are overlaid in Figure [| where it is seen 
that the structures and the absolute cross sections 
coincide. 



IV. ANALYSIS AND RESULTS 

The excitation function, shown in Figure [|, re- 
veals structure in the region from 1 MeV to 4 MeV 
that could be due to interfering broad resonances. 
A reasonable first assumption is that the struc- 
ture corresponds to the three lowest states in n N. 
This assumption is justified by the closed proton 
p 3 / 2 sub-shell in 10 C and agrees with the known 
predominantly single particle nature of the lowest 
states in 11 Be |20|], the mirror nucleus of 11 N. Tak- 
ing this as a starting point, we assume that the ob- 
served levels in n N are mainly of SP nature. The 
SP assumption validates the use of a shell-model 
potential to describe the experimental data of 11 N. 



A. Analysis of the three lowest levels in N 

The n N states are all in the continuum and the 
aim of the analysis was to obtain I* and other res- 
onance parameters as it can be done in the frame- 
work of the optical model. Because of the absence 
of other scattering processes than the elastic scat- 
tering channel, no imaginary part is included in 
the potential. The potential has a common form 
consisting of a Woods-Saxon central potential and 
a spin-orbit term with the form of a derivative of 
a Woods-Saxon potential. The Woods-Saxon (is) 
term has the usual parameters Vq (Vi s ), r (rg s ) 
and ao (ag s ) for well depth, radius and diffusity, 
respectively. Centrifugal and Coulomb terms were 
also included in the potential. The Coulomb term 
has the shape of a uniformly charged sphere with 
radius r c . The full potential is given in eq. [l^, 
where /i is the reduced mass of the system and A,r 
denotes the pion Compton wavelength. 
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(16) 



R = roA 1 / 3 , R c = r.A 1 / 3 , R ls = r A^ 3 

As a starting point, standard values of the po- 
tential parameters were chosen [5lJ and the well 
depths were varied separately for each partial wave 
(^=0,1,2), see Table |a. The cross-section of the 
experimental data at 180° was found to be larger 
than predicted by the potential model. As can be 
seen in the experimental spectrum, Figure ^, there 
is substantial amount of cross section above 4 MeV, 
indicating additional resonances in this energy re- 
gion. 
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FIG. 4. Experimental excitation function of n N to- 
gether with the best fit from the SP-model program. 
The energy is given as excitation energy of 11 N. The 
fit is made using a potential with parameters shown 
in Table |:c. The underestimation of the cross section 
around 1.8 MeV is the only failure of the potential 
model. This part is better described when the influ- 
ence of higher resonances are taken into account, see 
Figure []. 

The underestimation of the potential model can 
thus be attributed to influence of higher lying res- 
onances. To simulate the presence of those highly 
excited states, an amplitude / was added to the 
amplitude calculated from the potential model. 
The form of this extra amplitude was / = E(j b _ E , 



7 



where E$ was taken as a constant (4.5 MeV) and b 
was used as a parameter. As is seen in Figure |], the 
introduced amplitude is small in comparison with 
the measured cross sections, but it nonetheless was 
useful in the fitting procedure. A more sophisti- 
cated way to include the influence of higher reso- 
nances is to use an R-matrix procedure, and some 



attempts in this way were also made, see sec. IV B 



The best fit for conventional parameter values, 
only varying Vi is obtained for the level ordering 
lsi/2, Opi/2, and 0^5/2, and the corresponding pa- 
rameters are given in Table |:a. The curve result- 
ing from these parameters does not differ signifi- 
cantly from the one obtained using the parameter 
set Table |c. 

A potential with conventional parameters and 
the same well depth for all I will generate single 
particle levels in the order 0pi/2, lsi/2 and 
above the sub-shell. However, all attempts 

to describe the experimental data keeping this or- 
dering of the levels failed. A typical example of 
a calculated excitation function with this level se- 
quence is shown in Figure ||, with parameters in 
Table |b. This result is not surprising when con- 
sidering the well known level inversion in n Be. 




12 3 4 

Earn. (MeV) 

FIG. 5. A typical fit of the experimental data when 
the l/2~ state is assumed to be the ground state of 
11 N. It is evident that the excitation function is not well 
described. Putting the p state this low makes it too 
narrow, at the same time as the s state becomes very 
broad since it is now well above the Coulomb barrier. 

For the potential in Table |:c, the cross section 
of each partial wave is shown in Figure || together 
with the total calculated curve. Comparing the 
partial cross sections with the total cross section, 
it is clear that interference between the partial 
waves determine the shape of the total curve. The 
corresponding phase shifts are shown in Figure 0. 
The most common definition of resonance energy 
is where the phase shift 8e passes ir/2. As is seen in 
Figure 0, the phase of the l/2 + resonance, which is 
the broadest level, does not reach tt/2. Therefore, 



we have defined the resonance energy as where the 
partial-wave amplitude calculated at r = 1 fm has 
its maximum. The width is defined as the FWHM 
of the partial wave. One can note that for the 
l/2~ and 5/2+ levels, the same resonance ener- 
gies are obtained by our definition and 5i = tt/2. 
All attempts to change the resonance spins and 
parities or the order of the levels resulted in obvi- 
ous disagreement with the experimental data. We 
thus conclude that the unambiguous spin-parity 
assignments for the lowest states in 11 N are a 1/2+ 
ground state, a first excited l/2~ level and a 5/2 + 
second excited state. All further fitting procedures 
were performed with the aim to obtain more exact 
data on the positions and widths of the resonances. 



1000 




Earn. (MeV) 

FIG. 6. The partial waves s 1( / 2 , P1/2 



and d 



5/2 



to- 



gether with the total calculated excitation function for 
the best fit to the experimental data (Table |l|c). 
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FIG. 7. The phase shifts from the theoretically cal- 
culated curve which is the best fit to the experimental 
data (Table |c). 

A disadvantage of the potential model is that 
it produces resonances with single-particle widths. 
In general, the nature of the states is more compli- 
cated and their widths can be smaller than what 
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is predicted by the potential model. To investigate 
how changing the resonance widths would affect 
the overall fit, we changed the radius parameter ro, 
and fitted new well depths to get the best possible 
agreement with the data. It was evident that the 
widths obtained with ro = 1.4 fm are too large for 
the l/2~ and 5/2 + resonances, while r = 1.0 fm 
makes these levels too narrow. 

TABLE I. The sets of potential parameters used to 
obtain the best fits of the experimental data, and the 
resulting resonance parameters. Ve 3 =5.5 MeV is kept 
the same in all fits. The change of this parameter gave 
only minor modifications. The parameter b given in 
the footnotes is the value used in /(_B)=b/(4.5-E). 





Potential parameters 

Ve r , ri s , r c a a ls 
(MeV) (fm) (fm) (fm) 


Resonance 
(MeV) (MeV) 


° 1*1/2 

0f>3/2 
0^3/2 
0rf 5 /2 


-66.066 1.20 0.53 0.53 

-4z.ooD L.dU U.Oo U.OO 

-42.084 1.20 0.53 0.53 
-78.792 1.20 0.53 0.53 
-64.092 1.20 0.53 0.53 


1.30 1.24 
1.90 U.00 
-1.06 

4.40 0.90 
3.72 0.61 


° 1*1/2 
Ovi /o 

0^3/2 
0d 3 /2 

0d 5/2 


-45.360 1.40 0.65 0.30 
-33.474 1.40 0.28 0.30 
-32.340 1.40 0.53 0.30 
-58.086 1.40 0.53 0.30 
-45.570 1.40 0.35 0.30 


1.70 3.49 
1.11 0.11 
-1.22 

4.45 1.23 
3.75 0.60 


° 1*1/2 
Opi/2 
0P3/2 
0rf 3 /2 

0d 5/2 


-47.544 1.40 0.65 0.30 
-31.500 1.40 0.55 0.30 
-32.592 1.40 0.53 0.30 
-57.960 1.40 0.53 0.30 
-45.570 1.40 0.35 0.30 


1.27 1.44 
2.01 0.84 
-1.33 

4.5 1.27 
3.75 0.60 


a 1*1/2 
0pi/2 
0P3/2 
0rf 3 /2 
0rf 5 /2 


-56.280 1.20 0.75 0.30 
-42.420 1.20 0.55 0.30 
-42.210 1.20 0.53 0.30 
-78.960 1.20 0.53 0.30 
-62.874 1.20 0.50 0.30 


1.32 1.76 
2.14 0.88 
-1.33 

5.0 1.39 
3.79 0.59 


6 1*1/2 
Opi/2 
0P3/2 
0rf 3 /2 
0rf 5 /2 


-66.066 1.20 0.53 0.53 
-42.084 1.20 0.53 0.53 
-42.084 1.20 0.53 0.53 
-64.092 1.20 0.53 0.53 
-64.092 1.20 0.53 0.53 


1.30 1.24 
2.04 0.72 
-1.06 

9.87 4.53 
3.72 0.61 


1 1*1/2 
Ofl/2 
0P3/2 
0rf 3 /2 
0rf 5 /2 


-46.074 1.40 0.70 0.30 
-30.492 1.40 0.70 0.30 
-32.550 1.40 0.53 0.30 
-57.960 1.40 0.53 0.30 
-42.378 1.40 0.70 0.30 


1.27 1.56 
2.01 1.09 
-1.31 

4.50 1.27 
3.75 1.08 



a ro=1.2 fm, only varying Ve (6=1.25) 

b No ievel inversion (fe=1.25) 

c ro=1.4 fm, varying a and V( (b=1.25): 

the best fit to the data 
d ro=1.2 fm, varying a and Ve (b=2A) 
e The parameters used in |lj (6=0) 
f The parameters used to obtain the widths 

the singie particle limit. 



As the 1/2+ state is least affected by the change 
of radius the conclusion for this level is difficult, 
but the largest obtained width seemed most ap- 
propriate. Therefore, the radii parameters ro, rg s 
and r c were chosen as 1.4 fm, and the well depths 
and diffusities were varied separately for each £ to 
obtain the best fit of the experimental data up to 
4 MeV. An additional argument for choosing the 
larger radius was the fact that this parameter value 
gives a good simultaneous description of the mir- 
ror pair 11 Be and n N Q|. The curve obtained in 
this way that agreed best with the experimental 
excitation function is shown in Figure ^, and the 
corresponding potential parameters and resonance 
energies and widths are shown in Table |[c. For 
comparison, the values for ro=r^ s =r c =1.2 fm are 
also given in Table |:d. 

The extracted resonance parameters show a re- 
markable stability against changes in the potential 
parameter sets, meaning that different sets of pa- 
rameters that fit the data give similar resonance 
energies and widths. This is seen in Table |, com- 
paring different sets of parameters. The final en- 
ergies and widths are listed in table IV. The error 
bars include systematic errors and are dominated 
by a contribution from the spread in results from 
different parameter sets. Contributions from back- 
ground subtraction and solid angle corrections will 
be much smaller than those sources. 

The SP reduced widths could be extracted for 
the three lowest levels where the only possible de- 
cay channel is one-proton emission to the ground 
state of 10 C. The values of reduced widths are usu- 
ally presented as a ratio to the Wigner limit, which 
serves as a measure of the single particle width |l4| . 
In our case we have a way to give a more exact 
evaluation of the reduced widths as the ratio of 
the widths obtained from the shell-model poten- 
tial that fits the data (Table fic) to the widths cal- 
culated from a true shell-model potential. These 
ratios arc free from the uncertainties related with 
different definitions of the level widths. 

Since the true shell model potential is not known 
for n N, and we approximated this potential with 
the parameters shown in Table |:f. 

TABLE II. Reduced widths for the observed states 
obtained from the ratio of the widths from parame- 
ter sets c (experimental widths) and f (single particle 
widths) in Table [j| 



Level r 


Reduced width T exp /Fsp 


l/2 + 


0.92T0.2 


1/2- 


0.77T0.2 


5/2+ 


0.56T0.2 


3/2" 


0.15T0.2 
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Justification for using this particular set as shell 
model potential is that it simultaneously repro- 
duces the level positions in both n Be and n N 
and gives a width of the 1/2+ state that is larger 
than for the parameters in Table |:c. The reduced 
widths obtained in this way are given in Table g 

The potential parameters for the fit of the data 
at 180° (parameter set c in were used to describe 
the excitation function obtained by a detector at 
12.5° relative to the center of the chamber. The 
experimental data from this case are shown in Fig- 
ure H together with the theoretical curve without 
the scaling amplitude /. Comparing the exper- 
imental excitation functions in Figures |J and ||, 
rather big differences are seen. 
This reflects the fact that the laboratory angle de- 
pends on the interaction point in the chamber. The 
angular range goes from # c . m . = 150° for protons 
from higher resonances to 9 c . m =93° for low energy 
protons. This is taken into account in the calcula- 
tion of the excitation function, and from Figure || 
it is evident that the potential model describe the 
observed changes in the excitation function with 
angle, a fact which supports our interpretation. 



500 




-cm 



(MeV) 



FIG. 8. Experimental spectrum of protons in the 
detector placed at 12.5° relative to the center of the 
chamber. The full drawn curve is the result from 
the potential model using the parameters given in Ta- 
ble |:c. 



B. Analysis of the full excitation function 

In an attempt to investigate the influence of 
higher lying states on the cross section in the lower 
part of the experimental spectrum, a simplified R- 
matrix approach was used. For 11 Be, about 10 
levels are predicted in the energy region 2.7 MeV 
to 5.5 MeV |22] ], but only four resonances have 
been experimentally found |2C| ]. The knowledge 
of the levels in n N is even more incomplete, and 
our experimental data are not sufficient for a de- 



tailed R-matrix analysis. Therefore, the treatment 
described below was performed rather to outline 
possible questions than to give definite answers. 
The analysis was made using the potential model 
and adding resonances at energies above 4 MeV 
according to eq. [l7| 



35" = '"" ' = 



v-^E <' +1 > 1) c mtt+ " 

£ [e 21 ^ - l) e 2l (^+ CT «) 



(17) 

Two 

known levels in 11 Be (2.69 MeV, T = 200 keV, and 
3.41 MeV, r = 125 keV) were taken into account. 
The energy of those resonances in n N were deter- 
mined by calculating the Coulomb differences be- 
tween the mirror nuclei using the potential model. 
To fit the experimental data, the resonance ener- 
gies were varied around the value determined from 
the Coulomb-energy calculation. The values finally 
used in the R-matrix fit are shown in Table HI. 



±. 



Of. 



stands for states with j = £±1/2, resp. 
potential model amplitude at 180°, 
using the potential in eq. ^ 
resonance phase 
number of resonances 

phase relative to the hard sphere scattering 
for a particular resonance 
Coulomb phase of wave I 



The estimates of the widths of these states in n N 
are based on the known widths of the analog states 

TABLE III. The resonances used in the simplified 
R-matrix treatment. 

Potential model fit 





Potential parameters 


Resonance 




Vt r =res=r c 




E r 


R- 




(MeV) (fm) 


(fm) 


(MeV) (MeV) 


1/2+ 


-66.554 1.2 


0.5 


1.45 


1.56 


1/2- 


-41.286 1.2 


0.6 


2.13 


0.89 


5/2+ 


-64.801 1.205 


0.38 


3.74 


0.45 


Resonances added in the R-matrix fit 


*3/2 + 






3.94 


0.58 


3/2- 


Mirror of 2.69 MeV level n Be 


4.33 


0.27 


•(5/2+) 


Mirror of 3.41 MeV level 11 Be 


4.81 


0.40 


•(7/2-) 






5.4 


0.25 



The parameters for these states are only suggestions 
which reproduce the observed cross section. 
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in Be. Inclusion of these states already ac- 
counts for the missing cross section up to 3.7 MeV, 
but the part above 3.7 MeV is still underestimated. 

In particular, the energy region around 1.8 MeV 
is now better reproduced, indicating that interfer- 
ence of higher lying states indeed give the cross 
section that is not reproduced by the potential 
model in this region. Inclusion of a 3/2 + state at 
3.94 MeV and a high-spin state improves the de- 
scription also at energies above 4 MeV, as is seen 
in Figure ^. The parameters for the potential and 
included resonances are given in Table III. 

The conclusion drawn from comparing the re- 
sults in Table | and Table [II is that the posi- 
tions and widths obtained using only the poten- 
tial model are rather insensitive to the inclusion 
of higher states, which only modifies the absolute 
cross section. 



.a 
E. 

c; 




Earn. (MeV) 

FIG. 9. The fit of the potential model with added 
resonances. The full drawn curve includes the two 
known resonances in 11 Be plus a broad 3/2 + level 
around 4 MeV and a high-spin state at 5.4 MeV, which 
gives a good description of the full excitation function 
up to 4.7 MeV. Especially it should be noted that the 
theoretical curve now better reproduces the data at 
1.6-2.0 MeV. 

Of the three resonances included in the calcula- 
tions, only the one at 4.33 MeV is distinctly seen 
in the excitation function, see Figure ^. Its posi- 
tion corresponds to the 3/2 _ state at 2.69 MeV 
in 11 Be within 150 keV, and the cross section sup- 
ports a spin of 3/2 for this resonance. The obtained 
width of 270 keV also agrees with the width of the 
2.69 MeV level in 11 Be if decay by a l=\ proton is 
assumed. We thus conclude that the narrow res- 
onance at 4.33 MeV in n N is the mirror of the 
2.69 MeV state in n Be, both having I n =3/2~. 
The other resonances above 4 MeV are introduced 
in order to reproduce the cross section at higher en- 
ergies. The experimental data is not sufficient to 



give conclusive determination of any parameters of 
these states, but the existence of resonances above 
4.4 MeV is necessary to reproduce the measured 
cross section. 



V. DISCUSSION 
A. The three lowest resonances in N. 

Table || presents the parameters used in differ- 
ent fits of the deduced excitation function in the 
10 C+p system. The fits were all made under the 
assumption of three low- lying resonances. From 
these data we conclude that the three lowest states 
in n N have V= 1/2+, 1/2" and 5/2+. This is 
the first time all these states are identified in one 
single experiment. However, there have been in- 
dications of them in other reaction experiments. 
In the pioneering work on n N by Benenson et 
al. U), where the 14 N( 3 He, 6 He) reaction was stud- 
ied, it was proposed that the resonance observed 
at 2.24 MeV (r=740 keV) was a l/2~ state. This 
conclusion was based on the reaction mechanism 
in their experiment. Our data confirms this re- 
sult and both position and width are within the 
experimental errors of the two experiments. The 
difference may probably be attributed to different 
approaches in extracting the resonance parame- 
ters. In a recent paper by Lepine-Szily et al. |23|] , a 
state at 2.18 MeV was observed and interpreted as 
a l/2~ state, but with a width that was consider- 
ably narrower than in our work or that of Ref. |l3| . 

The state at 1.27 MeV, which we interpret as 
a l/2 + state, could not be seen in the two ex- 
periments in Refs. [^3|j2^], as the selected reac- 
tions quench the population of this state consid- 
erably. It could, however, be observed in an exper- 
iment performed at MSU where Azhari et al. [f24| 
studied proton emission from n N produced in a 
9 Be( 12 N, 11 N) reaction. They found indications of 
a double peak at low energies and by fixing the 
upper part of it to the parameters from Ref. []l3| 
they arrived at an excitation energy of 1.45 MeV. 

The 5/2+ (3.75 MeV) state was discussed in (§. 
The experiment presented in |p3| showed a state 
at 3.63 MeV with a width about 400 keV. The 
position of the resonance is close to ours but again 
the width is smaller in p3[ . 

As well as for the l/2 + and l/2~ states, the spin 
value for the 5/2 + resonance does not leave any 
doubt that it is the mirror state of the 5/2 + level 
at 1.78 MeV in n Be (T = 100 keV). The potential 
model with the parameters used for n N and given 
in Table |:c agrees very well for the width while 
the excitation energy becomes 1.63 MeV. Still we 
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consider this as an additional support for our in- 
terpretation. 

Fortune et al. pjJ have predicted the splitting 
between 0d 5 / 2 and lsi/2 states in U N from the sys- 
tematics of this energy difference for light nuclei, 
mainly assuming isobaric spin conservation. Their 
results can therefore be considered as a direct ex- 
trapolation of experimental data. The energy dif- 
ference obtained in our work (Ai?=2.48 MeV) is 
close to the prediction of 2.3 MeV in (^5j. The 
energy difference between the Is 1/2 and Opi/2 was 
calculated using the complex scaling method in ||(| 
and the value of 830 keV agrees with our data 
which gives 740 keV. 



B. Resonances above 4 MeV 

We interpret the structure around 4.3 MeV as 
due to a sharp resonance in n N, which is the mir- 
ror state of the 2.69 MeV level in n Be, Table fll. 
Several different experiments (see for example |22fl ) 
give the spin-parity for this 11 Be level as 3/2~. 



The negative parity is well established from mea- 
surements of the 11 Li beta-decay [p7|-p9[, and by 
measurements of the 9 Be(t,p) 11 Be reaction p^ |. 
There is also good agreement between the Cohen- 
Kurath prediction for the spectroscopic factor and 
the reduced single particle widths of these mirror 
states. We found very good agreement between the 
widths if the states undergo nucleon decay with 
£ = 1 (J = 3/2). If the states decay with orbital 
momentum, 1=2 (j = 3/2), the state in n N will 
be at least twice as broad, and in the case of £ — 3 
(j = 5/2) it would be at least 3 times broader. Also, 
for £ = 3 the reduced single particle width will be 
too large, contradicting pPf] . Using the simplified 
R- matrix approach, the position of the 3/2 _ level 
was determined as 4.33 MeV. The observed cross 
section for the population of this state is also in 
accordance with a 3/2" assignment. The calcula- 
tions further indicate that about one third of the 
width of the 4.33 MeV state is due to the pro- 
ton decay to the first excited state in 10 C. Even a 
small branch of this decay results in a large reduced 
width. This indicates a large coupling of 3/2~ 
state to the first excited state of the core, as was re- 
cently predicted by Descouvemont plj . In Q we 
proposed a different structure for the 3/2 _ state 
(two particles in the Is state) because preliminary 
treatment resulted in a too small width (70 keV) 
for the state. 

In the present experiment there is an experimen- 
tal cutoff at 5.4 MeV (see Figure |^) and the exci- 
tation function increases towards this high-energy 
end. This is most likely due to higher-lying states 



but we cannot make any assignments for them 
based on our data. However the authors |24|] had to 
introduce a broad (r= 500-1000 keV) state in the 
energy region around 4.6 MeV to explain the spec- 
trum from 11 N decay. They proposed the broad 
state to be a 3/2~ state. Our data show that the 
3/2~ state in n N is rather narrow, and therefore 
another state has to be assumed in order to explain 
the data in [ p4j . This is also a justification for the 
inclusion of the 3/2 + resonance is our R-matrix fit. 

Various theoretical calculations (for recent ref- 
erences see J||]) have attempted to reproduce the 
level sequence in 11 Be. Most models emphasize the 
role of coupling between the valence neutron and 
the first excited 2 + state in 10 Be in generating the 
parity inversion. It is well known that model as- 
sumptions influence the wave functions more than 
their energy eigenvalues and thus models giving the 
correct level sequence predict very different core 
excitation admixtures. For the ground state in 
11 Be, the admixtures given by theoretical calcu- 
lations vary from 7% §l||§ to 75% @. Theo- 
retical results are frequently compared to spectro- 
scopic factors obtained from nucleon transfer re- 
actions. The single-particle spectroscopic factors 
for 11 Be have been obtained from 10 Be(d,p) reac- 
tions fjjj. Even if the theory of stripping reac- 
tions is very well developed, many parameters are 
involved in the extraction of these results from the 
data. Evaluating single-particle nucleon widths us- 
ing a potential model involves fewer parameters. 
For the lowest states of n N we obtained the re- 
duced widths given in Table [nj. For the S!/ 2 state 
we have a reduced width of ~1 which, taking the 
15% experimental error in the width into account, 
indicates that no large core-excitation admixtures 
are needed to describe the ground states of n N 
and n Be. 



VI. SUMMARY 

The excitation function in the 10 C+p system 
has been studied using elastic resonance scattering. 
The low-energy part was analyzed in a potential 
model while the high-energy part was described in 
a simplified R-matrix approach. The ground state 
and the first two excited states in the unbound 
nucleus n N was found to have the spin-parity se- 
quence of l/2 + , l/2~ and 5/2 + which is identical 
to that found in its mirror partner n Be. A narrow 
3/2~ state at 4.33 MeV was identified as the mirror 
state of the 2.69 MeV state in 11 Be. The energies 
and widths of the observed states are listed in Ta- 



ble [V . The agreement among experiments as well 
as between our results and theoretical calculations 
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TABLE IV. A summary of all experimental and theoretical results on n N. E r and F r denote the resonance 
energy and the width of the resonance, respectively. The error bars in this work include the systematic errors (25 
keV in the cm. frame) as well as the spread in results obtained for different potential parameters that all fit the 
data. 





1/2+ 






1/2- 




5/2+ 




Ref. 


E r 








r r 







Experimental papers 



-I 07 +U.l« 

-0.05 



-0.1b 

-0.05 



This work 



1.44 ±0.2 



2.01 



0.84 ±0.2 



3.75 ±0.05 



0.60 ±0.05 



[23 
[ 24[ a 

Theoretical papers 



1.45 ± 0.40 



> 0.4 



2.24 ±0.1 
2.18 ±0.05 
2.24 ±0.1 



0.74 ±0.1 
0.44 ± 0.08 
0.74 ±0.1 



3.63 ±0.05 



0.40 ± 0.08 



[35 
[25 
[36 

[37 



1.54 
1.60 ± 0.22 
1.4 
1.1 
1.2 



0.62 

, ,+1.0 

1.31 

0.9 
1.2 



2.24 ±0.1 
2.49 
2.21 
1.6 
2.1 



0.74 ±0.1 
1.45 
0.91 
0.3 
1.0 



3.74 
3.90 
3.88 
3.8 
3.7 



0.3 
0.88 
0.72 
0.6 
1.0 



a For the 1/2 state the results from 
'The results obtained with tq = 1.45 fm is presented. 



are very satisfactory. 

The quasi-stationary character of n N states was 
used to evaluate the reduced single-particle widths 
for the identified states. This result indicates 
small coupling between the valence nucleon in the 
ground state n N and the first excited 2 + state in 
10 C, and the same conclusion should be valid for 
n Be. 

The experimental technique to use elastic- 
resonance scattering with radioactive beams has 
proven to be a very efficient tool for investigations 
beyond the driplinc. 
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